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ABSTRACT 


A review of the governing equations and boundary con- 
ditions used to model the response of ablative materi- 
als submitted to a high-enthalpy flow is proposed. The 
heritage of model-development efforts undertaken in the 
1960s is extremely clear: the bases of the models used in 
the community are mathematically equivalent. Most of 
the material-response codes implement a single model in 
which the equation parameters may be modified to model 
different materials or conditions. The level of fidelity 
of the models implemented in design tools only slightly 
varies. Research and development codes are generally 
more advanced but often not as robust. The capabilities 
of each of these codes are summarized in a color-coded 
table along with research and development efforts cur- 
rently in progress. 
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NOMENCLATURE 


Latin 

T{ Diffusion flux of the i th species, kg • m~ 2 • s ~ 1 
m Mass flow rate, kg • m~ 2 • s -1 

Ai Gaseous species i 

Aj Arrhenius law pre-exponential factor, SI 

Ch Stanton number for heat transfer 
Cm Stanton number for mass transfer 
c p Specific heat, J • kg~ x • K~ x 
e Specific energy, J • kg~ 2 

Ej Arrhenius law activation energy, J - kg~ 2 
Fj Fraction of mass lost through pyrolysis reaction 
3 

h Specific enthalpy, J - kg 1 
j Diffusive flux, mol • m~ 2 • s _1 

Ki Chemical equilibrium constant for reaction i 

l Thickness or length, m 

mj Arrhenius law parameter 

Molar mass of species k, kg • mol* 


N g Number of gaseous species 

rij Arrhenius law parameter 

N p Number of pyrolysis reactions 

q Heat flux, J • m~ 2 • s -1 

R Perfect gas constant, J • kg~ x • AT -1 

v Convection velocity, m • s -1 

y Mass fraction 

Fo Forchheimer number 

K Permeability 

p Pressure, Pa 

Greek 

[3 Klinkenberg coefficient, Pa 

e Volume fraction 

7 ji Stoichiometric coefficient, reaction j species 

g Viscosity, Pa • s 1 

(j Reaction rate, mol • m~ 3 • s _1 

uj s Solid reaction rate, mol • m~ 3 • s _1 

n Pyrolysis gas production rate, kg • m~ 3 • s -1 
p Density, kg • m ~ 3 

r Characteristic time, s 

fj Advancement of pyrolysis reaction j 

Subscripts 

a Ablative material (gas, fiber, and matrix) 
c Char 

e Boundary layer edge properties 

/ Reinforcement (non-pyrolyzing phase) 

g Gas phase 

m, PM Polymer matrix 
mv Virgin polymer matrix 

p Pyrolysis 

pg Pyrolysis gas 

s Solid phase 

Conventions 

d x -() Divergence 

dt() Time derivative 

T Second order tensor 

u Vector 



1. INTRODUCTION 


During re-entry, a fraction of the heat is transferred to 
the thermal protection system (TPS) leading to a grad- 
ual temperature increase of the material (figure 1). With 
the temperature increase, the virgin material is succes- 
sively transformed and removed by two phenomena. The 
first transformation phenomenon is called pyrolysis. Dur- 
ing pyrolysis, the pyrolyzing phase of the material (of- 
ten a polymer matrix) progressively carbonizes and loses 
mass producing pyrolysis gases. The pyrolysis gases are 
transported out of the material by diffusion and convec- 
tion through the pore network. During this transfer, their 
chemical composition evolves as their temperature in- 
creases. The second transformation phenomenon is the 
ablation of the char that is composed of the residual car- 
bonized matrix and of the non-pyrolyzing phase (often 
a carbon or silicon-carbide fibrous preform). Depending 
on reentry conditions, ablation may be due to heteroge- 
neous chemical reactions (oxidation, nitridation), phase 
change (sublimation), and/or mechanical erosion (spalla- 
tion). Material response models should predict accurately 
the ablation rate and the peak temperature of the bondline 
at the interface of the TPS and the substructure. The first 
open literature publication providing a very detailed and 
comprehensive analysis of ablative material response in 
high enthalpy environments is the Aerotherm report from 
1968 describing their suite of design tools [16]. CM A 
and ACE, Aerotherm’s in-depth material response and 
surface ablation codes respectively, are cited as a refer- 
ence in most publications in the field. The models imple- 
mented in current design-capable tools are mostly repli- 
cas (or parallel developments) of the Aerotherm model 
with slight variations. Interestingly, the modifications to 
the Aerotherm model have mainly involved simplifica- 
tions, with some of the Aerotherm capabilities currently 
no longer maintained in most of the major tools. Recent 
interest in manned-rated and challenging design missions 
(e.g. high mass, very high velocity, porous materials) has 
raised the need for high-fidelity models capable of pro- 
viding optimized design and comprehensive uncertainty 
quantifications. All the capabilities of the Aerotherm’s 
suite of tools and the rich academic work on pyrolysis 
and ablation are being revisited and progressively intro- 
duced (or re-introduced) both in research codes and in 
design tools. In a complementary effort, several aca- 
demic, government, and industrial teams are working on 
the development, the implementation, and the validation 
of original physics-based models that will enable anchor- 
ing of CMA/ACE-based design tools, accurate uncer- 
tainty analysis, and maybe become the future base mod- 
els for design-rated codes. This paper first presents the 
pyrolysis-ablation problem through the governing equa- 
tions (mass, momentum, and energy conservation) and 
boundary conditions. Different levels of modeling fi- 
delity are presented and discussed. An effort was made 
to gather information on the simulation tools that are ac- 
tively used either for design or for research and develop- 
ment. The capabilities of each of these codes are summed 
up in figure 4. 
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Figure 1. Picture of a core of ablative material extracted 
from the TPS of Stardust [33] and schematic of the zones 
of degradation illustrating the material response to a 
high-enthalpy flow. 


2. CONSTITUTIVE MODELS 


2.1. Mass conservation 

The gaseous mass-conservation equation includes a pro- 
duction term (right-hand side) to account for the pyrolysis 
gas production, noted II, and reads 


9t(e g p g ) + d x -(e gPg v s ) = II (1) 


In several codes, the time derivative is omitted and the 
gas flow problem is treated as a succession of steady state 
problems (see section 4). This simplification is correct 
when the variation of the intensive variables (tempera- 
ture, pressure) are slow compared to the characteristic 
time of the flow in the porous medium. The character- 
istic time of the pyrolysis gas flow, r pg may be defined 
as the ratio of the thickness of the char layer - l c - to 
the velocity of the gas. In typical re-entry applications, 
r p g = l c /vg ~ 0.01/1 = 0.01s. Therefore, the omission 
of the time derivative is an acceptable practice for situ- 
ations for which the variations of the intensive variables 
are negligible over time steps of r pg T ste p = Is. The 
determination of the direction of the gas velocity, v g , is 
necessary to solve the average mass-conservation equa- 
tion. In several one-dimensional codes, this equation is 
numerically integrated with the assumption that the gas 
flow is perpendicular to the surface and directed towards 
the surface. This is exact in one-dimensional steady-state 
problems with an impermeable back face. In other con- 
ditions and in multi-dimensional problems, the direction 
of the flow has to be determined by resolution of the 
momentum-conservation equation. 

The pyrolysis gas production is obtained by fitting ther- 
mogravimetry analysis of the resin decomposition using 



one or several Arrhenius laws. For example, for pheno- 
lic polymers, it has been shown that the pyrolysis degra- 
dation process follows four steps [36], that may be de- 
scribed by four heterogeneous decomposition reactions 
[37]. A convenient notation for j G [1, N p \ pyrolysis re- 
actions is 


iv 9 

PMj -»• T: 7 jjAj (2) 

i= 1 

where PMj is a Active solid species of the pyrolysing 
polymer matrix (PM). The pyrolysing matrix density is 
then given by 


This requires the experimental determination of the stoe- 
chiometric factors - 7 ji, which are not directly available 
in the literature but may be derived from experimental 
studies [34, 36, 37]. The overall pyrolysis gas production 
may still by obtained from: II = J2i=li [^iMi\. 

Higher fidelity models are being developed and imple- 
mented. They account for species production, trans- 
port, and chemical reactions (finite-rate chemistry) within 
porous media. The species conservation equation may be 
written in mass fraction - yi as 


^ti^gPgVi) + (tg PgViV g) — 7T Mi EgidiMj, 

( 8 ) 


^rnPrn 


N v 


€ mvPmv 




3 = 1 


(3) 


where 


(1 - £)”*' 


T nj Aj exp 



(4) 


The pyrolysis-gas production is given by 


N p 

n = ^ti^mPm) = e mvPmv ^ ^ Pj &t (Jjj ) (3) 

3 = 1 

In the literature, the equations used to describe pyrolysis 
models vary but they are mathematically equivalent. 

It is important to mention that state-of-the-art design 
codes do not track the species production. Only the aver- 
age mass production ~n - is computed from the Arrhenius 
laws. A constant elemental fraction of the pyrolysis gas 
is assumed. The gas chemical composition and derived 
quantities (gas enthalpy, viscosity, mean molar mass) are 
then computed using the chemical equilibrium assump- 
tion or heuristic methods. 

The pyrolysis gas production rate for each species i could 
readily be obtained using 


N p 

(^"i) ^ mPmv Y [dtO F i 7 ji] (6) 

3 = 1 

where 


7 ji 


7 a 


2/c=l IjkMik 


(7) 


Both pyrolysis species production - 77 - and chemical 
species production - uj t - are needed. For the compu- 
tation of uoi, the finite-rate chemistry model developed 
by Pike and April in the late 1960s [28, 2] is used for 
preliminary analyses. This model was developed using 
chemical data and experimental techniques available at 
this time. Efforts are being undertaken in several teams to 
develop up-to-date finite-rate chemistry models and mod- 
em experimental set-ups for their validation [ 20 ]. Ti is 
the diffusion flux of the i th species. At low pressures, 
mass transfer (diffusion) in porous media is not negligi- 
ble compared to advection [18]. Mass transfer in porous 
media is a complex problem. The effective diffusion co- 
efficient is smaller than the bulk diffusion coefficient due 
to the porosity and the tortuosity of the material [18]. 
Several models are available. A popular extension of the 
Stefan-Maxwell model [12] to porous media is the dusty 
gas model [25]. To our knowledge no fully capable ab- 
lation material-response code has such a capability yet. 
High-fidelity models including in-depth finite-rate chem- 
istry will need to account for diffusion. 

The solid-phase mass conservation is also integrated to 
compute the effective density of the solid. The volume- 
averaged density change of the matrix (due to pyrolysis 
-n) is currently modeled using forms equivalent to 


dt{e m Pm ) = n (9) 

which is easily derived from equation 5. The current as- 
sumption is that there is no ablation or coking in-depth. 
Coking is neglected and ablation is modeled as a surface 
phenomenon only. Therefore, ablation is accounted for 
using a prescribed recession velocity at the wall, handled 
as a boundary condition (rather than as an in-depth con- 
stitutive equation) as described in section 3. 

Current research efforts aim at developing models for 
in-depth coking and ablation. For this application, the 
solid mass-conservation equation may be generalized to 
account for heterogeneous reactions 

d t (e s p s ) = d t (e m p m + e f p f ) = -II ~Y e g UJ i M i 


( 10 ) 



However, the determination of the heterogeneous reac- 
tion rates -w£- is not an easy task and is still being inves- 
tigated. 


2.2. Momentum conservation in porous media 

The average gas velocity is obtained by resolution of the 
momentum-conservation equation. In porous media, the 
volume- averaged momentum conservation may be writ- 
ten as 


v 


g 


l \ - Vp 

Cgfi 1 + Fo 


K'd t p 


(ID 


Most of the materials are anisotropic, therefore, the per- 
meability - K - is a second order tensor. For exam- 
ple, Fiberform, the carbon preform of PICA [35], has 
orthotropic permeability properties [23]. For creeping 
(Stokes) flows in the continuum regime (in the pores of 
the material), the volume-average momentum conserva- 
tion degenerates into Darcy’s law (/? = 0, Fo = 0). The 
term 1 + /3/p is the Klinkenberg correction to account for 
slip effects (at the pore scale) when the Knudsen number 
(ratio of the mean free path to the mean pore diameter) is 
not small. The term 1 + Fo is the Forchheimer correction 
to account for high velocity effects at the pore scale (flow 
separation in the continuum regime). Typically, Forch- 
heimer effects are expected to occur for pyrolysis gas ve- 
locities higher than 50 m/s (that is, in high-density ab- 
lative materials submitted to very high heat fluxes). It 
is not advised to use both corrections simultaneously as 
they address different regimes. 


2.3. Energy conservation 


According to Puiroux [29], solid and gas phases are 
in thermal equilibrium as long as the Peclet number 
for diffusion of heat within the pores is small (Pe = 
e gPg c p,gd P v g/k g ). In most of the applications of interest 
for space agencies, the small pore size (< 100 pm) and 
the slow pyrolysis gas flow (v g ^ 1 m/s) insure a small 
Peclet number: the gas temperature accommodates to the 
solid temperature within the pores. Under the thermal 
equilibrium assumption, the energy conservation may be 
written as 


Pa&a — ^gPg^g "F ^ mPmhm ~F ^fPfhf (13) 


The second and third terms of the left-hand side are the 
energy convected (advection) and the energy transferred 
(diffusion) by the pyrolysis gases, respectively. Heat 
transfer is conveniently modeled as an effective diffu- 
sive transfer (Fourier’s law). The effective conductivity 
- k - is a second order tensor accounting for conduction 
in the solid, conduction in the gas, and effective radia- 
tive heat transfer. The validity of this volume- averaged 
approach is questionable. The main issue is the valid- 
ity of the linearization of the radiative heat transfer. A 
theoretical study has shown that radiative heat transfer 
may be linearized for two-dimensional carbon-fiber pre- 
forms [40, 41]. The applicability to other materials is not 
straightforward and needs to be investigated. The second 
term on the right-hand side is the energy dissipated by 
viscous effects in Darcian regime [10]. It is small com- 
pared to the heat transfer term and often neglected. 


3. BOUNDARY CONDITIONS 


At the bondline, conservative boundary conditions are 
generally used (adiabatic and impermeable). At the wall 
and in ablative conditions, surface energy balance and 
surface mass balance are used as boundary conditions. 
[Of course, simple wall boundary conditions may always 
be used for simple analyses, e.g. fixed temperature.] 


3.1. Surface energy balance 
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Figure 2. Energy balance at the wall 


The surface energy balance at the wall depicted in figure 
2 reads 


Qconv (pV)h w "F Fad, in Fad, out Fond 


~\~tfipghpg H - rn ca h ca — 0 


(14) 


Ng 

dtPa^a "F d-x.’i^gPghgV g) "F <9 X ' 

i=t 

= 5x-(k-5xT) + ^(K _1 -v)-v 


where the total (storage) energy of the ablative material 
is the sum of the energy of its components 


where the convective heat flux - q con v = PeU e C' H (h e — 
h w ) - and the radiative heat flux are extracted from CFD 
simulations. The Stanton number Cr is corrected to ac- 
count for the blockage induced by the pyrolysis-ablation 
gas-blowing; that is, the heat transfer coefficient is cor- 
rected. For example, the following correction is widely 
used C' H = C#/n(l + 2XB')/ln(2XB '), where B' = 
(' rhpg + rn ca ) / (p e u e C m) is a dimensionless mass flow 
rate and A is a scaling factor usually taken equal to 0.5 
[6]. 


3.2. Surface mass balance and recession rate 


The pyrolysis-gas flow rate - m pg - is directly obtained 
in the material-response code by integration of the py- 
rolysis, transport, and mass equations, as explained pre- 
viously. However, the ablation rate - rh ca - is a func- 
tion of both the mass transfer in the boundary layer and 
the thermo-chemical properties at the wall (pyrolysis- 
gas blowing rate and composition, temperature, pressure, 
boundary-layer gas composition). A common practice 
is to assume thermochemical equilibrium at the wall to 
compute the ablation rate. The model still in use in the 
community was developed in the sixties [27] . It is based 
on element conservation in steady- state in a control vol- 
ume close to the wall as sketched in figure 3. The un- 
derlying hypothesis is that over a time increment At, the 
equilibrium chemistry problem in the control volume is 
quasi- steady (decoupling of the material response and of 
the boundary layer problem). This increment At should 
be at least as long as the time increment of the heat trans- 
fer simulation (material response code) but short enough 
so that p, T, rh P g , and y pg variations may be neglected. 
This is verified in typical applications. For this presenta- 
tion, we shall assume equal diffusion coefficients of the 
elements. Failure modes (spallation, mechanical erosion) 
are not included and the char is assumed to be composed 
of a single element (for example, carbon). 
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Figure 3. Element mass-fraction conservation at the wall 


The inputs and outputs to this problem are: 


• Since p, T are fixed, the element mass-fraction con- 
servation in the control volume is equivalent to the 
mass conservation. 


Under the hypotheses that Prandtl = Lewis = 1 and 
that the diffusion coefficients are equal for the elements, 
equation 15 may be rewritten as 

Pe^eC H (Vk,w yk,e)3~(pV)yk,w = ™pgyk,pgP'dlcayk,ca 

(16) 

where, Ch is the Stanton number and (pV) = m pg + 

rh ca . 

The formation reaction of species Ai may be written: 

A% x ^ ^ V i,kAfc (17) 

kEElements 

The i chemical equilibriums read: 

y Vi,kln(x k ) - ln(xi) - ln(Ki) = 0 (18) 

kEElements 

with xi — 1 if Ai is a solid species. Species mole frac- 
tions sum to one: 


y Xi = i (i9) 

iESpecies 


To sum up, the set of equations solved is: 

Pe'U'eC H (]Jk,w yk,e)F^pV^)yk,w — Al pg y k ^ pg FTTl ca y k ^ ca 

( 20 ) 

y Vi,kln(x k ) - ln(xi) - ln(Ki) = 0 (21) 

kEElements 

with X{ = 1 if Ai is a solid species. 

E Xi = 1 (22) 

iESpecies 


• Inputs: rh pg , y k , pg , y k , ca = 1, y k , e , p , T. 

• Outputs: m ca , y k , w . 

The conservation of the mass-fraction of element k in the 
control volume close the the wall reads: 

jk,w (pV)y k , w = kh pg y k , pg H- rri ca y k , ca (15) 


The base model may be extended when needed to ac- 
count for multicomponent mass transfer, non-equal dif- 
fusion coefficients, failure (spallation, melting), a solid 
phase made of more than one element (example: SiO 2 ), 
corrections to account for heterogenous finite-rate chem- 
istry. Current development efforts aim at fully modeling 
the boundary layer and coupling it to material codes, with 
the recession directly computed in the flow solver. 

The surface pressure is an input to the material code (ob- 
tained from CFD simulations). It is the boundary condi- 
tion required for the averaged momentum equation. 


where pg= pyrolysis gases, ca = char ablation products, 
w= wall (or control volume). Notes: 


4. SIMULATION TOOLS 


• The relative mass fractions sum to 1 in each phase A , ^ u w . ^ ^ , 

r More than twenty ablative-material response tools are 

J2k yk,w = 1; J2k yk^pg = 1; Y.k yk^ca = 1 currently in use or in development for hypersonic re-entry 


applications. The name of the codes and useful infor- 
mation are provided in table 1 . The contact listed is ei- 
ther the code developer or a current active user. For each 
code, one open-literature reference is provided (for most 
of them, many are however available). Our understand- 
ing of the current code capabilities and/or development 
strategies based on open-literature publications are sum- 
marized in figure 4. 


5. CONCLUSION 


A review of the governing equations and boundary con- 
ditions used to model the response of ablative materi- 
als submitted to a high-enthalpy flow has been proposed. 
The heritage of model-development efforts undertaken 
in the 1960s is extremely clear: the bases of the mod- 
els used in the community are mathematically equivalent. 
Most of the design-rated material-response codes imple- 
ment a single model in which the equation parameters 
may be modified to model different materials or condi- 
tions. The level of fidelity of the models implemented in 
design tools only slightly varies. Research and develop- 
ment codes developed for analysis - at least in a first stage 
- are generally more advanced but often not fully capable. 
To sum- up, more than twenty codes are currently in use 
or in development, with an active community both main- 
taining state-of-the-art capability and seeking to increase 
the fidelity of the state-of-the-art model. 
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Table 1. List of currently available simulation tools 


Name 

Contact 

Owner 

Users 

Applications 

References 

Amaryllis 

T. van Eekelen 

Samtech, Belguim 

EADS Astrium, ESA 

Design 

[39] 

CAMAC 

W.-S. Lin 

CSIST, Taiwan 

Taiwan Ins. of Sci. Tech. 

Unknown 

[21] 

CAT 

N. N. Mansour 

NASA ARC, USA 

NASA ARC 

Analysis 

[22] 

CHALEUR 

B. Blackwell 

SNL, USA 

SNL 

Design 

[5] 

CHAP 

P. Keller 

Boeing, USA 

Boeing 

Design 

[14] 

CMA 

R. Beck 

Aerotherm, USA 

NASA, SNL 

Design 

[15] 

CMA/SCMA 

C. Park 

Tokyo Univ., Japan 

JAXA 

Design 

[1] 

CMA/KCMA 

P. Reygnier 

ISA, France 

ISA/ESA 

Analysis 

[30] 

COYOTTE2 

D. W. Kuntz 

SNL, USA 

SNL 

Design 

[17] 

FABL 

J. Merrifield 

Fluid Grav. Eng. Ltd., UK 

ISA/ESA/FGE 

Analysis 

[32] 

FIAT 

Y.-K. Chen 

NASA ARC, USA 

NASA, SpaceX 

Design 

[6] 

FIAT3D 

Y.-K. Chen 

NASA ARC, USA 

NASA ARC 

Analysis 

[7] 

HERO 

M. E. Ewing 

ATK, USA 

ATK 

Analysis 

[11] 

ITARC 

M. E. Ewing 

ATK, USA 

ATK 

Design 

[11] 

libAblation 

R. R. Upadhyay 

Univ. of Tex. Aust., USA 

UTA 

Analysis 

[38] 

MIG 

S. Roy 

Univ. of Flo., USA 

Univ. of Florida 

Analysis 

[4] 

MOPAR 

A. Martin 

Univ. of Mich., USA 

UKY/Univ. of Michigan 

Analysis 

[24] 

NEQAP 

J. B. Scoggins 

N. Carol. St. Univ., USA 

NCSU 

Analysis 

[31] 

NIDA 

G. C. Cheng 

Univ. Alab. Birm., USA 

UAB 

Analysis 

[42] 

PATO 

J. Lachaud 

NASA ARC, USA 

Univ. Calif. Santa Cruz 

Analysis 

[19] 

PRESENT 

J. Dec 

NASA LaRC, USA 

NASA LaRC 

Analysis 

[9] 

STAB 

B. Remark 

NASA JSC, USA 

Fluid Gr. Eng. 

Design 

[8] 

TITAN 

F. S. Milos 

NASA ARC, USA 

NASA 

Analysis 

[26] 

TMU 

A. R. Bahramian 

T. Modares Univ., Iran 

TMU 

Analysis 

[3] 

US3D 

G. Candler 

Univ. of Minn., USA 

UM 

Analysis 
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Gas-phase M ass Conservation 

In-depth : Eq. 1 

Storage (a t ...) 
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Pyrolysis model 

In-depth: Eq. 2-7 



SoA Arrhenius laws (->n) 
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Species production (->ni) 

























Gas-species Conservation 

In-depth: Eq. 8 

Storage (at ...) 
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M ulti-component diffusion (dxF) 
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Finite-rate chemistry ( ni , cos) 
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Solid-phase mass conservation 

In-depth: Eq. 9-10 

Pyrolyzing matrix mass loss 


























In-depth ablation/coking 
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Momentum conservation 

In-depth: Eq. 11 

Darcy's law 


























Klinkenberg/Forchheimer 
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Energy conservation 

In-depth: Eq. 12-13 

Storage (at ...) 
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Effective conduction 


























Viscous dissipation 


























Boundary conditions 

At the wall: Eq. 14-22 

Surface energy balance 
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Wall chemistry from B' table 
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Internal wall chemistry solver 
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Other utilities 

| Integrated libraries 

Equilibrium chemistry solver 
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Integrated boundary layer code 
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Script -coupling to CFD code 
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Figure 4. Simulation- tool list and capabilities. [The authors wish to apology in advance for any missing or incorrect 
information contained in this figure. Corrections and addenda will be greatly appreciated .] 
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. Introduction 

Physics and Chemistry in Ablative Materials 



Surface 

Char 

layer 


Pyrolysis 

zone 


Virgin 

material 

i 

i 

i 


Type 1 model 



Q] 


Core - Stardust TPS [1] 


[1] M. Stackpoole et al., Post-Flight Evaluation of Stardust Sample 
Return Capsule Forebody Heatshield Material, AIAA 2008-1202 


Color code for the presentation 

Three types of material-response 
models have been identified & 
defined [during the 4th AFOSR / 
SNL/N ASA Ablation Workshop. 
March 1-3, 2011, Albuquerque, 
NM]: 

• Type 1 : Heat transfer, pyrolysis, 
simplified transport of the pyrolysis 
gases, equilibrium chemistry, 
surface ablation (current state-of- 
the-art) 

• Type 2: Type 1 + averaged 
momentum equation for the 
transport of the pyrolysis gases 
(e.g. Darcy’s law) 

• Type 3: High-fidelity model 


TType 2 + Animate <^wM^ raphs m 
multi-component diffusion, in-depth 
ablation/coking, explicit radiative 
transfer model, ...) 



. Outline 

Motivation: Pure curiosity, review currently available codes and models 


1 - In-depth modeling 

1.1 - Mass conservation 

1.2 - Mass transport 

1.3 - Energy conservation 

2 - Boundary conditions 

2.1 - Surface energy balance 

2.2 - Surface mass balance and recession rate 

3 - Summary 

3.1 - Table 1 : List of codes, contacts, and references 

3.2 - Table 2: Codes capabilities 
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1.1. Mass conservation (gas/solid system) 

Candid illustration of the chemistry of Pyrolysis & Ablation : incomplete mechanism! 


Shock 


20 



Ablation 

(leads to surface 
recession) 

Chemistry of the 
pyrolysis gases 


Pyrolysis & 

Pyrolysis-gas 

production 
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1.1. Mass conservation (gas/solid system) 

Models: Type 1 / Type 2 / Type 3 - In green: Experiments - 

Type 1 codes - hypothesis: instantaneous transfer (no storage) 


<v(/&)= n 

i_ 


Pyrolysis gas flux 


Pyrolysis laws (Arrhenius) = f(Temperature 
- Thermogravimetry analysis (TGA) 


Type 2 & Type 3 

a *(vO +3 *'(v>A)- n 


Open porosity 
- Pycnometry 



see 

Momentum 

Conservation 


v 

Mean molar mass:Af = \ X f M i 


Need to 
determine Gas 
Composition 
(in mole 
fractions of the 
species) 



Equilibrium chemistry 
+ elemental 


composition 


Finite-rate chemistry 
+ species conservation 


1.1. Species conservation (gas/solid system) 

Models: Type 1 / Type 2 / Type 3 - In green: Experiments 


• Gaseous Species Conservation 

d t ( E gP g y t )+ d x • ( £ g p g y, v , )= + £ g ^,M i 


i 


I 


Mass fraction of 
species i 


See mass 
transport 


I 

Finite-rate chemistry mechanism 
- Flow reactor + spectro/chromato 


Pyrolysis production rate for each species 
-Thermogravimetry analysis (TGA) + mass spectroscopy / chromatography 


• Solid-phase volume-fraction conservation (in-depth ablation/ 
coking) 

Porosity increase/decrease due to heterogeneous reactions 


■ hypotheses: riq swelling, constant intrinsic solid densities.-- 

Fibers: o f e , = -S2, > sco. Carbonized matrix : d,£ = -12 

t f ^ f / j i t m m 


t f ^ f 

Molar volume 


£CO ; 


Quantification of porosity change/heterogeneous chemistry 
-Scanning Electron Spectroscopy (SEM) / Tomography before and after tes 


ting. 



1.1. Species conservation (gas/solid system) 

Pyrolysis-gas SPECIES-production modeling 




1.1. Species conservation (gas/solid system) 

Finite-rate chemistry of the pyrolysis gases 



[April69] April, G. C., “Energy Transfer in the char Zone of a Charring Ablator”, 

Ph. D. Dissertation, Department of Chemical engineering, Louisiana State University, 1969. 
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1.1. Mass conservation (gas/solid system) 



p = 0.12 atm ; T=898 K -> Recession: 0.7 cm in 1 hour. 

Two principal observations: 

- reduction of fiber diameter eventually leading to recession 

- sharp ablation front (about 0.2 mm) [i.e. oxygen is quickly consumed, reaction » transpor 


Validation of a volume-averaged fiber-scale model for the oxidation of a carbon-fiber preform. 

J Lachaud, N N Mansour, Joe Marschall, et al. To appear, 42nd AIAA Thermophysics Conference, 27-30 June 2011. 


1.2. Mass transport in porous media 

Models: Type 1 / Type 2 / Type 3 - In green: Experiments - 


Species velocity 



Average gas velocity 
Momentum Conservation in 
porous media 


Diffusion velocity 

Multicomponent mass transfer in porous media 

(e.g. Dusty Gas Model: Stefan-Maxwell + Bosanquet) 


r 1 K 1 + fi/p 

v = — d x p 

g £ n 1 + Fo 

o 

Generalization of Darcy’s law 


K: permeability 

(3: Klinkenberg coefficient (low 
density, 

slip effect) -typically p < 0.1 atm 
Fo: Forchheimer number (high 
velocity, 

separation) - typically v g > 50 m/s 
_ Permeameter 

from the 9 a 

composition, cf options for M 


s 


In porous media, the effective diffusion 
coefficient is reduced (because wall 
collisions decrease the effective mean free 
path) 




D 


ref 




Tortuosity 

- Diffusivity apparatus 

- Tomography + DNS 
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1.3. Energy conservation 
Models: Type 1 / Type 2 / Type 3 



- In green: Experiments - 


Thermal equilibrium between the phases studied by N. Puiroux [AIAA 
2002-3336]. His conclusion is that the thermal equilibrium assumption 
is valid. 


One temperature model 


PA = e g p g e + s 



Computed or 
neglected 


Solid enthalpies 
-Differential Scanning 
Calorimetry (DSC) 


Effective conductivity 
(conduction + radiative transfeif) 

- Guarded Hot Plate 

- Flash diffusivity 
-Arc Jet calibration 


ii 
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2.1. Surface Energy Balance 

At the wall 



Convective flux 

1 

Advective flux 

A 

Radiative heating 
| 

Radiative cooling 

1 

ablating 

^conv 

T (PV)H W 

1 qradin 

1 ^radout - 

surface 


1 ^cond 

1 m pg H pg 

t m H 

ca ca 


Conduction flux Pyrolysis gas (pg) flux 

Char ablation (ca) flux 
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2.2. Surface mass balance and recession rate 

At the wall (simplified description: equilibrium chemistry, no spallation) 


Wall 


Mass transfer flux 


Advection flux 



ca 


Pyrolysis gas (pg) flux Char ablation (ca) flux 


1- Conservation of element mass fraction in the control volume (y k w ) 

jk,w “1“ ( )l/k,w — rilpgl/k,pg H - ^ caVk,ca 


2- Equilibrium chemistry 



Provides both the ablation rate (m ca ) and the gas composition (y k w ) 
in the control volume (and derived quantities, like wall enthalpy) 



The dimensionless tabulation of m ca = f(m pg , T, p) is called a B’ tables 


. Outline 

Motivation: Pure curiosity, review currently available codes and models 


1 - In-depth modeling 

1.1 - Mass conservation 
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1.3 - Energy conservation 
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2.1 - Surface energy balance 

2.2 - Surface mass balance and recession rate 
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3.1 - Table 1 : List of codes, contacts, and references 

3.2 - Table 2: Codes capabilities 
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3.1. Table 1: List of codes, contacts, and references 

We found 25 codes actively used and with open-literature references... 


Name 

Contact 

Owner 

Users 

Applications 

References 

Amaryllis 

T. van Eekelen 

Samtech. Belguim 

EADS Astrium. ESA 

Design 

[17] 

CAMAC 

W.-S. Lin 

CS1ST. Taiwan 

Taiwan Ins. of Sci. Tech. 

Unknown 

[18] 

CAT 

N. N. Mansonr 

NASA ARC. USA 

NASA ARC 

Analysis 

[19] 

CHALEUR 

B. Blackwell 

SNL, USA 

SNL 

Design 

[20] 

CHAP 

P. Keller 

Boeing. USA 

Boeing 

Design 

[21] 

CM A 

R. Beck 

Aerotherm. USA 

NASA, SNL 

Design 

[22] 

CMA/SCMA 

C. Park 

Tokyo Univ.. Japan 

JAXA 

Design 

[23] 

CMA/KCMA 

P. Reygnier 

ISA. France 

ISA/ESA 

Analysis 

[24] 

coyotte: 

D. W. Kuntz 

SNL. USA 

SNL 

Design 

[25] 

FABL 

J. Merrifield 

Fluid Grav. Eng. Ltd., UK 

ISA/ESA/FGE 

Analysis 

[26] 

RAT 

Y.-K. Chen 

NASA ARC. USA 

NASA. SpaceX 

Design 

[27] 

FI ADD 

Y.-K. Chen 

NASA ARC. USA 

NASA ARC 

Analysis 

[28] 

HERO 

M. E. Ewing 

ATK. USA 

ATK 

Analysis 

[29] 

ITARC 

M. E. Ewing 

ATK, USA 

ATK 

Design 

[29] 

libAblation 

R. R. Upadhvay 

Univ. of Tex. Aust.. USA 

UTA 

Analysis 

[30] 

MIG 

S. Roy 

Univ. of Flo.. USA 

Univ. of Florida 

Analysis 

[31] 

MOPAR 

A. Martin 

Univ. of Mich.. USA 

UKY/Univ. of Michigan 

Analysis 

132] 

NEQAP 

i. B. Scoggins 

N. Carol. St. Univ.. USA 

NCSU 

Analysis 

[33] 

NIDA 

G. C. Cheng 

Univ. Alab. Birm., USA 

UAB 

Analysis 

[34] 

PATO 

J. Lachaud 

NASA ARC. USA 

Univ. Calif. Santa Cruz 

Analysis 

[35] 

PRESENT 

J. Dec 

NASA LaRC. USA 

NASA LaRC 

Analysis 

[36] 

STAB 

J. Merrifield 

NASA JSC, USA 

Fluid Gr. Eng. 

Design 

[37] 

TITAN 

F. S. Milos 

NASA ARC. USA 

NASA 

Analysis 

[38] 

TMU 

A. R. Bahramian 

T. Modares Univ.. Iran 

TMU 

Analysis 

[39] 

US3D 

G. Candler 

Univ. of Minn.. USA 

UM 

Analysis 

[40] 


... please help us extending the list! 
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Code capabilities 
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Solid-phase mass conservation 

In-depth: Eq. 9*10 

Pyrolyzmg matrix mass loss 
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Momentum conservation 
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Energy conservation 

In-depth: Eq. 12-13 
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Boundary conditions 

At the wall: Eq. 14-22 

Surface energy balance 
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Integrated boundary layer code 



P 

















p 






Script -coupling to CFD code 
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Table 2 : 

Codes capabilities 
Please accept our 
apologies for any 
error. 

Corrections and 
suggestions of 
improvement will 
be greatly 
appreciated! 

Conclusion 

- 25+ codes actively 
used 

- same base model 

(from the 1960s) 

- design-rated tools: 
simple and robust 

- analysis and 
research 
tools: generally 
more 

advanced but often 
not fully capable for 
design 

- very active and 
diverse community 
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1.3. Energy conservation 

Radiative heat transfer in fibrous preforms. 2D DNS & volume averaging 


left wall 


periodicity 


right wall 




Question : is it correct to linearize the 
radiative heat transfer in fibrous 
media? 

Method : DNS using geometrical optics 
[Refraction and diffraction negligible for 
T > 700K (i.e. when radiation not 
negligible)] 

Expression of the effective 
onductivi4v fr aT 3 Ax: 


- Collision based ray-tracing algorithm 
compute the form factors 

- Finite-element code SAMCEF 


corn 

u,m 


Ifo^jfig Ljne^rizatjpnxiWhe radiative 
transport) a ' extgr 


Objective : compute the extinction 
coefficient and verify the accuracy of the 

Numerical Validation of an Effective Radiation Heat Transfer 

T. van Eekelen, J. Lachaud. Journal of Spacecraft and Rockets. Technical note. Accepted for 
publication. 
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T (K) 


2.3. Sensitivity analysis 


2011 ablation-modeling test-case (4th AFOSR/NASA/SNL Ablation workshop, March 1-3, 2011) 
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14 participants : we are still working on the analysis of the results, but a few remarks 

• 13 codes (type 1 and/or type 2) gave results in agreement with FIAT (cf. PAM1, PAM 2, 
CAT). Some differences attributed to: numerical errors (+/- 1%) and misinterpretation of the data 
for 2 codes (+/-6%). 

• 1 code out-of-range (difference > 200%) 

• Type 3 (with finite-rate chemistry) predicts lower temperatures (PAM 3). 


Mass fractions, Yi (-) 


2.3. Sensitivity analysis 

2011 ablation-modeling test-case (4th AFOSR/NASA/SNL Ablation workshop, March 1-3, 2011) 


Ablation Workshop Test-case 2011 - Pyrolysis-gas composition 
Code: PATO/PAM_3 (type 3 - finite-rate chemistry) 



surface (x=0) 


x (m) 
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